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PREFACE 


A photochemical model for the Venus atmosphere including 
oxygen, hydrogen, chlorine and sulfur chemistry is presented. 
Primary emphasis is placed on the formulation of the sulfur 
chemistry and its role in the formation of H 2 S 0 ^ molecules 
for the hydrated sulfuric acid droplet cloud deck of Venus. 

The present study both complements and extends the 
earlier sulfur chemistry of Prinn (1973, 1975) and the hydrogen, 
chlorine and oxygen chemistry of Sze and McElroy (1975) . 

Sulfur is released by photolysis of COS and oxygen is supplied 
by photolysis of CO 2 . Subsequent reactions oxidize sulfur 
and lead to formation of H 2 S 0 ^. 

The dominant sulfur oxidation process for H2S0^ production 
involves reactions of S with O 2 followed by reactions of SO and 
SO 2 with atomic oxygen. The odd hydrogen production scheme 
proposed by Prinn (1973) plays a relatively minor role. The 
abundance of O 2 controls the rate limiting step in the H 2 SO^ 
formation sequence. The concentration of SO 2 and COS calculated 
from the model are comparable above 65 km. The upward flux of 
COS is balanced by the downward fluxes of SO 2 and H2S0^ in 
approximately equal proportions. The supply of H2S0^ may be 
sufficient to maintain a 10 km thick cloud deck if the sulfuric 
acid cloud droplets have a residence time of about a few years. 
Shorter residence time requires either a thinner cloud or 


(iii) 


P|i£CC0U^ p®; 


additional chemical reactions to increase the efficiency of 
converting SO2 into SO^. 

Further studies of the chemical model are warranted to 
explore the sulfur chemistry. Future measurements on SO2 
abundance would be extremely useful in refining the model 
results . 
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CHAPTER I 


■% 


INTRODUCTION 


Understanding of the composition and structure of the atmos- 
phere of Venus has progressed significantly in recent years. Both 
earth based astronomy and spacecraft measurements have made impor- 
tant contributions. These measurements have provided observational 
constraints for development of atmospheric models. Photochemical 
reaction schemes used in these models have to date placed primary 
emphasis upon the stability of the CO 2 atmosphere. Several mech- 
anisms, involving catalytic chemistry of odd hydrogen (H, OH, 

HO 2 , H 2 O 2 ) and odd chlorine (Cl, CI 2 , HCl, CIO, ClOO) , were suc- 
cessfully employed by McElroy, Sze and Yung (1973) and Sze and 
McElroy (1975) , for recombination of CO 2 in the Venus atmosphere 
above the visible H2S0^ cloud layer. Alternative mechanisms 
have been proposed by Prinn (1973, 1975) , who stressed sulfur 
photochemistry for production of H 2 S 0 ^ molecules above the cloud, 
and recombination of CO 2 below the cloud. A more complete photo- 
chemical model for the Venus atmosphere would combine the merits 
of the oxygen, hydrogen and chlorine chemistry of Sze and McElroy 
with the sulfur chemistry of Prinn. Development of such a numer- 
ical model and initial evaluation of it t're the subjects of this 
report. 

Interest in exploring this more complete photochemistry 
for Venus has increased with the realization that updates chem- 
ical reaction rates have lowered the efficiency of the catalytic 
recombination schemes presented by Sze and McElroy (1975) . With- 
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out other oxygen sinks such as that provided by sulfur, the excess 
0 , produced would violate the observed mixing ratio, unless it 
could be transported downward more rapidly and converted to CO2 
below the cloud layer. Efforts to specify the fraction of the 
excess ©2 that can be transported below the clouds and the frac- 
tion that can be converted into H2S0^ above the clouds, are 
reported. 
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CHAPTER II 

PHOTOCHEMISTRY FOR VENUS 
2.1 Structure and Composition of the Atmosphere 

The vertical structure and composition of the Venus atmos- 
phere determined by earth-based and spacecraft measurements is 
rather poorly defined. The major constituent of the atmosphere 
is CO 2 with a coliunn abundance equal to or in excess of 93% 

( Vinogradov et al., 1971). Nitrogen and argon are likely 
candidates for the other major gases that make up most of the 
remaining 7%. The surface pressure and temperature of Venus are 
approximately 93 bars and 750° respectively ( Marov et al., 1973). 
The surface is however completely obscured from view by a cloud 
layer composed of partially hydrated sulfuric acid droplets 
which is located approximately in the 45-62 altitude range. 
Vertical profiles for the number density of CO<, and temperature 
from the cloud tops to a 100 km altitude are given in Table 1. 

The mixing ratios of trace gases which have been detected 
in the Venus atmosphere or which have measured upper limits, are 
listed in Table 2. The small concentration of O 2 (<lppm) and CO 
('SOppm) were initially most surprising results, given the ease of 
producing large amounts of CO and ©2 from photolysis of CO 2 . This 
sparked interest in understanding the stability of the CO 2 atmos- 
phere. The trace species H 2 O is much more concentrated below the 
Venus cloud layer, suggesting that much of the upward flow of 
wate,r is both photochemically converted into H 2 SO^ and absorbed 
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Table 1 


Number Density and Temperature of the Venus Atmosphere 


Height NCsl 

(km) (cm"^) 


T(z) 

(°K) 


62 

6.80(18) 

255 

64 

4.24 (18) 

246 

66 

3.17(18) 

239 

68 

2.10(18) 

232 

70 

1.55(18) 

225 

72 

9.99(17) 

217 

74 

7.26 (17) 

209 

76 

4.53(17) 

200 

78 

3.22 (17) 

194 

80 

1.90(17) 

188 

82 

1.30 (17) 

184 

84 

8.17 (16) 

180 

86 

4.50(16) 

177 

88 

2.90(16) 

175 

90 

1.73(16) 

173 

92 

9.91(15) 

171 

94 

6.30 (15) 

170 

96 

3.67 (15) 

168 

98 

2.00(15) 

167 

100 

1.27(15) 

166 


NASA SP-8011 Models of Venus Atmosphere (1972) 


4 


Tabxts 2 


Abundance of Trace Gases in the Atmosphere 


Trace Gas 
(above clouds) 

Mixing Ratio 
(ppm) 

Reference 

CO 

50 

young (1972) 


<1 

Traub & Carleton (1974) 

H 2 O (above cloud) 

0.5-40 

Barker (1975) 

(below cloud) 

0 

I 

0 

Rossow & Sagan (1975 ) , 
Vinogrador et al. (1971) 

HCl 

0 

• 

CTV 

Young vl972) 

COS 

<1 

Cruikshank (1967) 

so^ 

<0-035 

Cruikshank & Kuiper (1967) 

H^SO^ (in clouds) 

*a,o 

Rossow & Sagan (1975) 

H 2 S 

^0 . 3 

Anderson et al. (1969) 

HP 

-0 . 01 

Young (1972) 

He 

10 

Kumar & Bx'oadfoot (1975) 
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by the cloud droplets. HCl through i.hotoIysis, providts a major 
sourc® of hydri-^iTtn and chlorine to the Venus atmosphere. This large 
concentration of HCl, together with the small concentration of 02» 
is responsible for the lack of an ozone layer in the Venus atmo'* 
sphere. The abundance of H 2 S 0 ^ molecules above the cloud layer is 
not known but the appro-ximate amount present below the cloud tops i® 
given. Upper limits are given for other sulfur species that are 
likely to be present in the atmosphere. Measured values are also 
given for KF and He, although neither is included in chemical 
modeling. 

The vertical structure and composition of trace gases in the 
atmosphere reflects an interplay between diffusion, cliemical pro- 
duction and chemical loss processes of the various species. Tur- 
bulent mixing, v;hich is modeled by eddy diffusion, dominates the 
vertical transport from the surface of Venus to the homopause 
(turbopause) , located at an altitude of about 130 km. The indivi- 
dual gas species in this region are driven to a constant mixing 
ratio in the absence of chemical production and loss. The role of 
global circulation in vertical transport is not clear but it may 
also be important below the homopause. Molecular diffusion dominates 
vertical transport from the homopause to the exosphere (at about 
200 km) , with ballistic transport operating above this altitude 
level. The lighter gas species therefore have greater relative 
abundances as the altitude increases, with chemical processes 
altering this behavior locally. 

Chemical modeling to be considered will focus attention upon 
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th® altitud® range from the cloud tops to 100 km. This is the 
region where sulfur is actively converted into H2SO^ molecules 
and transported downward to supply the cloud. Reaction schemes for 
oxygen, hydrogen, chlorine and sulfur chemistry to be used in 
the Venus atmosphere model are presented below. 

2.2 Pre-Sulfur Chemistry 

The hydrogen and chlorine chemistry for the Venus atmosphere 

was developed out of the need to understand the stability of the 

CO2 atmosphere. Figure 1 shows the photodissociation rate for 

CO2, HCl and H2O in the Venus atmosphere as calculated by McElroy, 

Sze and Yung ( 1973 ) . The column dissociation rate for CO5 is 
13 -2 -1 

about 1x10 molecules cm sec , a value capable of depleting 
Venus of its entire CO2 atmosphere in a few million years. 

Photochemical models using simple oxygen chemistry given 
by the reaction scheme (2.2.1) were unsuccessful in recombining CO^. 


CO^ ii * 00 + 0 


(2.2. la) 

0 + 0 + COt - 0., 

+ CQ^ , 

(2.2. lb) 

CO + JSO2 CO2 


(2.2. Ic) 

CO + 0 + M CO2 

+ M 

(2.2. Id) 


Reaction 

reaction 


(2.2. Ic) 
(2.2. Id) 


is endothermic at Venus temperatures, while 
is operative, but has an extremely small rate 
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constant because it is spin forbidden. 


A more complex chemistry involving other atoms was therefore 
considered to explain the stability of CO,. Successful catalytic 
recombination schemes for this purpose using hydrogen and chlorine 
chemistry were given by Sze and McElroy (1975). In their model, 
removal of CO proceeded by two main paths: 

H + O2 + CO2 ^ no 2 + CO2 

0 + HO2 " OH + O2 

CO + OH CO2 + H 

and 

2 t H + O2 + CO2 ^ IIO2 t CO2 > 

HO^ + HO2 ^ H2O, + O2 
hv + H2O2 ^ 20H 
2 1 CO + OH “ CO^ + H } 

The first corresponds to the net reaction 

H + CO + 0 H + CO, (2.2.4) 

while the second is equivalent to 

2H + 2CO + 0, > 2H + 2CO2 (2.2.5) 

The necessary source of hydrogen is provided by photolysis of 
HCl 

hv + HCl ^ H + Cl (2.2.6) 

followed by 

Cl + FI2 HCl + II (2.2.7) 


( 2 . 2 . 2 ) 
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ware molecular hydrogen consumed in (2.2.7) is produced above the 
cloud by 

H -»■ HO 2 Hg + O 2 (2.2. C) 

and supplied from below by an upward flux. Both rocombination 
schemes (2.2.2) and (2.2.3) depend upon the odd hydrogen species HO 2 . 
Similiar reaction schemes exist for chlorine using the radical 
CIO, but they are of lesser importance. A list of relevant reac- 
tions for oxygen, hydrogen and chlorine chemistry is given in 
Table 3. 

New kinetic data on chemical reactions available since 1975 
have shown that the rate of the reaction 

HO 2 + Cl HCl + O 2 (2.2.9) 

is about If.H) larger than the value adopted in the work of 

Ste and ''iCBlroy . The increase in this reaction rate reduces 

the efficiency of the catalytic recombination scheme in two ways. 
First it reduces the concentrations of HO 2 and Cl. Second it 
removes two chemically active species and replaces them with the 
more stable molecules HCl and 0-. The concentrations of 0, 0., and CO 
must then increase. To avoid violating observational constraints, 
the excess CO and 0, must be transported downward and converted 
below the cloud deck by unspecified thermochemical processes into 
CO.,. The CO, below the clouds is then recycled by upward trans- 
port. If however sulfur chemistry is added to the model, the 
excess O 2 may alternatively be able to locally supply oxygen 
for formation of H 2 S 0 ^ molecules. 
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2.3 Sulfur Chemistrv 


Formation of H 2 SO^ molecules in the Venus atmosphere may be 
limited by the supply of hydrogen, oxygen and sulfur atoms. In 
the sulfur chemistry to be considered, hydrogen is supplied above 
the clouds by photolysis of HCl and from below the clouds by an 
upward flux of H 2 and H 2 O. Oxygen is supplied by photolysis of 
CO 2 so that the downward flux of H 2 S 0 ^ must be less than the col- 
imm dissociation rate o£ COj 

r '^CO **00 “ 1.5x10^^ om"^ seo'^ (2.3.1) 

B 2 2 

Sulfur is supplied above the clouds by photolysis of COS, a post- 
ulated but undetected molecule in the Venus atmosphere. COS is 
postulated to be reformed by thermochemical processes below the 
cloud deck and returned by an upward flux. The upper limit on 
the mixing ratio of COS above the cloud layer is about lo””. 


Sulfur Reaction Scheme 

An explicit photochemical reaction scheme for production 
of H 2 S 0 ^ based on these ideas was first introduced by Prinn 
(1973). The Prinn scheme consists of six reactions; 


COS 

+ 

hv 

CO 

+ 

S('D) 

[ (a) 

S('D) + 

CO2 

s 



(b) 

S 

4* 

©2 

so 

+ 

0 

(c) 

SO 

+ 

O H 

so., + 

H 

(d) 

SO^ 

•f 

HO2- 

SO3 

H 

(e) 

S 03 


H 2 O 

+ CO2 

H2S0^ + CO2 

(f) 
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The rate limiting step in scheme (2.3.2) is the reaction 
(2.3.2e). The number density of HO2 is about 10 cm at 
cloud tops, dropping rapidly with increasing altitudes ac- 
cording to Sze and McElroy (1975) . The rate constant for 

—16 3 —1 

this reaction is very slow, with a value 9 x 10 cm s 
The conversion of SO2 to H2SO^ is therefore rather inefficient. 
A modest production rate for H2SO^ is insured only by requiring 
an appreciable concentration of SO2 which might however violate 
the observed upper limits of SO2 column abundance. 

An alternate sulfur reaction scheme may be summarized as 
follows. 


COS + hv 

-> 

CO + S('D) 

(a) 

S('D) + CO 2 


S + CO 2 

(b) 

S + 02 

- 

SO + 0 

(c) 

SO + 0 + M 


SO 2 + M 

(d) 

SO 2 + 0 + M 

->• 

SO 3 + M 

(e) 

3 t H 2 O + M 


H,SO^ + M 

(f) 


The above scheme involves only sulfur and oxygen species, in 

contrast to reaction scheme (2.3.2) which involves odd hydrogen, 

in addition to oxygen and sulfur chemistry. The rate limiting 

step in scheme (2.3.3) is the reaction (2.3.3e). The rate 

—33 6 —1 

constant for this reaction is about 7 x 10 cm sec , which 

-14 3 -1 

xs equivalent to a two body rate of about 5 x 10 cm sec , 
a moderately fast rate. 

The concentration of atomic oxygen is significantly 
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enhanced because of the major revisions in rate constants as 
discussed in Section 2.2. Moreover, the concentration of 0 
would be further enhanced by including the reaction 

S + O 2 -»■ SO + 0 (2.3.4) 

which represents an additional source of oxygen below 70 km. 
Figure 2 shows the importance of atomic oxygen production due 
to teaction (2.3.4) relative to CO 2 photolysis. It appears 
that scheme (2.3.3) might provide an interesting alternative 
to scheme (2.3.2), as proposed by Prinn (1973, 1975). 


The complete chemical model to be explored in this study 
is given in Table 4 and incorporates all reactions from both 
of these schemes. Additional sulfur reactions are also in- 
cluded. Rate constant subscripts are numbered consecutively 
with those in Table 3. The rate k 2 g is unmeasured and an 
estimated value is given. Reaction k^j^ is an effective re- 
action, a composite of the reaction 


SO 2 + OH + M -V HSO 3 + M 


(2.3.5) 


followed by photolysis of HSO^. The reaction Jg is not 
effective in the lower portion of the atmosphere where sulfur 
has its impact because CO 2 provides an optically thick path 
for photons below wavelengths of about 2100 S. Rate constants 
kg 2 kgg atc treated as parameters and simulate postulated 
heterogeneous production processes for H 2 SO^ . 
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Altitude (km) 


CO 2 hv CO+0 



S +O 2 SO + 0 


] 0 ^ 

Reaction Rate (cm"*' sec"') 


Figure 2 . 


Production Rate for Atomic Oxygen. The rate of 
creation of oxygen atoms -by photolysis of CO2 
and the reaction of S and ©2 are indicated as a 
function of altitude. 
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CHAPTER III 


MODEL DEVELOPMENT 


3.1 Atmospheric Model 

A one dimensional photochemical model similiar to the one used by 
fse and McElroy (1975), has been developed for the Venus' atmosphere. 
The number density for the ith species, n^^, is governed by the steady 
state continuity and diffusion equations; 


d-J 

dz 


i . p . - L.n. 


(3. 1.1. a) 


and 


^"i . "i + ;^i 


-K ( — + 
dz H 


av 


T dz 


"i dT 

-D. (— i + -i. + (1 +a.) -i — ) , O.l.l.b) 
dz T dz 

where P. and L. n. represent the production and loss rate (molecule cm”^ 

•U JL 

s”^) ; K and D^^ represent the eddy and molecular diffusion coefficients, 
and represent the average scale height and species scale height 
and is the effective thermal diffusion coefficient. 


The eddy diffusion coefficient was assiuned to have the form 

* 

Z < Z < Z 




(3.1.2) 


= K 


< Z 1 Z^ 
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where i«{Z) is the number density of the CO2 background gas and 

K* . K(Z*). Ths molecular diffusion coefficient Oj has the form 


0 ^( 2 ) 


M 


3 1 

nW) 


f 2kT(Z) .^ 

' u 1 ^ 

i 


(3.1.3) 


where is the rigid sphere cross section for collision of the 
tih 

i species with the background gas, and is the reduced mass 

t*]h 

for the i species mass m^^ and the background gas mass m. 

Vertical profiles for the number density n(Z) and temperature 
T(Z) of CO2 are given in Table 1 for altitudes from the cloud tops 
(62 km) to 100 km. 


The altitude dependence for the scale heights are given by 


H 


av 


Hi 


- ^^(2) 
rug 

- kT(Z) 


mig 


(3.1.4) 

(3.1.5) 


where g is the acceleration of gravity at the Venus surface. The 


coefficient is effectively zero except for hyUorgen. 


Eddy diffusion dominates molecular diffusion below the homo- 
pause (-v 130 km) , while molecular diffusion is dominant above 
this level to the exosphere. The lower boundary, Z^, is chosen 
to be at 62 km and the upper boundary Z^ at 100 km, below which 
sulfur chemistry is incorporated in the model. 
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3 . 2 Nuraarical Procadura 


Solutions for tha diffusion equation can be divided into 
two categories. 

(1) those for species involved in chemical reactions that 
are fast compared to diffusion processes and 

(2) those for species involved in chemical loss times 
which are only comparable or slower than the charac- 
teristic diffusion times. 

For the first category, the number densities may be obtained 
by solving a set of non-linear coupled algebraic equations. 

P • ** L j n . 

A generalized Newton ' s method was employed to slove these equations . 

In the later category the second order differential equation 
(3.1.1) for n^ must be solved. A finite difference method is 
applied and an iterative procedure adopted to obtain a solution 
for these species. Two boundary conditions are required for 
each equation and must be chosen in a consistent way so as not 
to overspecify or underspecify the solution for the set of equations. 

For the photochemical reaction schemes adopted for oxygen, 
hydrogen, chlorine (see Table 3) and sulfur (see Table 4) , there 
are 11 species in the first category and 7 species in the second 
category. The fast species are 0, 0^, S, SO, SO^, OK, H, H 2 O 2 , 

IIO 2 , CIO and Cl 2 . The slow species are O 2 , CO, HCl, Cl, COS, SO,, 
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I 


and Ii2S0^. The two remaining trace species H 2 and 112*^ assumed 

—7 

to hava constant mixing ratios of 2 x 10 and 10 respectively. 

The chemical loss terms are given in Table 5 for the photochemical 
reaction scheme of Table 3 and Table 4. 

3.3 Boundary Conditions 

The number density of eac»i slow trace gas species (see section 
-.3) is governed by the diffusion equation (3.1.1). For each species, 
two boundary conditions are required for a unique solution. Either 
the flux or mixing ratio of each species may be specified at the 
lower boundary (Z « Zg) . The flux is usually specified at the top 
boundary (2 =* 2,^,) . 

Boundary conditions must be chosen sc as to provide number 
density profiles consistent with observations. For Venus (see Table 
2) average column mixing ratios are known for CO and HGl. Oniy upper 
limit values have been measured for O 2 , COS and SO 2 . Boundary con- 
ditions are therefore not well defined for most of the slow species. 
More measurements are clearly needed to eliminate this difficulty. 

Solutions to be explored here are based upon the assumption 
that photolysis of CO 2 and COS above the cloud layer provide the 
major sources of oxygen and sulfur for H 2 SO^ formation, excluding 
an upward flux of H 2 O. It then follows that fluxes of CO, S02?02 and 
H2S0^ must be downward. Conservation of atoms then requiring an 
upward flux for COS, CO 2 and H20- 'Jhe fluxes of O 2 / CO, HCl, Cl, 

COS and SO 2 are taken to be zero at the top boundary (2 = 2^) . 


21 


Production and Loss Terms for Trace Species 
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Boundary conditions for the lower boundary are summarized in 


Table 6 , where f. and >J). denote the mixing ratio and flux of the 
ith species respectively. At this boundary, a range of values will 
be considered for the O 2 and COS mixing ratios. The mixing ratios 
of CO and HCl are chosen so as to agree with observations. The 
chlorine number density is assumed to have its chemical 

equilibrium value. The downward flux of SO^ is set equal to a frac- 
tion t of the upward COS flux, where 


0 < a < 1 


The diffusion equation (3.1.1) is not solved for the species 
H2S0^. To determine a number density profile for H2S0^ requires 
either specification of its mixing ratio at boundary (an answer to 
be hopefully deduced) or specification of a chemical loss process 
such as would occur by condensation of H 2 S 0 ^ above the cloud tops. 
The downward flux of H2S0^ may be calculated without specification 
of boundary conditions and is sufficient to evaluate the physical 
parameters to be defined in section 3.4. 

3.4 Discussion of the Physical Processes 

The major questions raised by the addition of sulfur chemistry 

to the oxygen, hydrogen and chlorine chemistry for Venus can be 

understood by considering the fate of oxygen atoms produced by 

photolysis of CO 2 above the Venus cloud layer. The column production 

2 

rate of oxygen (the number of oxygen atoms produced per cm per sec) 
above the cloud layer is given by 

Po = "coj 
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Boundax-y Conditions 
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Some of these oxygen atoms cycle through the hydrogen and chlorine 
catalytic recombination schemes and react with CO to reform C02* 
The column loss of 0 by these processes is given by 


^0 * ^002*^^ • (3.4.2) 

If the CO 2 recombination schemes is incomplete (i.e. Pq^Lq) ' 
the excess oxygen must be transported below the cloud tops by 
downward fluxes of oxygen bearing species, where, through un- 
identified thermochemical processes, it is converted into CO 2 and 
COS and recycled upward. In the chemistry under consideration, this 
excess oxygen may be removed by downward fluxes of 02 » SO 2 and H 2 S 0 ^. 
In addition to this an upward flux of H 2 O at the cloud tops may 
provide a partial source of oxygen in the production of 
through the reaction 


SO, + HoO + M '>■ H,S0^ + M 


(3.4.3) 


A conservation of oxygen atoms above the cloud tops then 
requires that 

Po ■ Lo + =0 (3.4.4a) 


where the fluxes and are directed downwards and 

O 2 SO 2 "2°°4 


the flux Q directed upwards. Dividing this expression by P^, 
we have 


^ “ ^cat ” ^02 ” ^S02 ’ ^H2S0^ ° 


(3.4.4b) 
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where 


X , * 
cat 

Lo/Pq 

(3.4.5a) 


(24«02)/Po 

(3.4.5b) 

’‘sOj = 


(3.4.5c) 



(3.4.5d) 


The quantity is the fraction of the oxygen, produced by photoly- 

sis of CO 2 above the cloud tops, which is catalytically reconverted 
to CO 2 above the cloud tops. For 100% conversion The quan- 

tities Xq , XgQ and Xjj gQ are the fractions of the oxygen produced 
2 2 2 4 

by photolysis of CO 2 which are transported below the cloud tops by 
downward fluxes of 02 » SO 2 and H 2 SO^ respectively. 


The efforts of Sze and McElroy (1975) , which included no sulfur 
chemistry, emphasized the case where and Xq -vO whereas the 

efforts of Prinn (1973,1975), with no hydrogen and chlorine chem- 
istry, emphasized the case where X^j g^ -)-l. Values for the parameters 

2 4 

^cat' ' ^SO so include the coupling between oxygen, 

hydrogen, chlorine and sulfur photochemistry are sought for the 
model given in Tabel 3 and Table 4 . 


The efficiency of the sulfur chemistry in producing H 2 SO^ 
molecules may be better understood by considering the fate of a 
MOfnv a+-nm nrnduoiad by nhohol'^sis of COS abo”"- 

This is accomplished by considering the general conservation relations 
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given in Table 7 for the 5 different atoms, 0, F, C, Cl and H 
which make uo all the chemical soecies under consideration. 

The quantity in Table 7 is the escape flux cl K at the Venus 
exosohere. Notice that the relation (3.4.4) derived earlier is 
nothing more than the conservation relation for oxygen atoms in 
Table 7, where and ‘I’^Qg have been eliminated bv the carbon 
conservation relation. 


To explore the sulfur budget, we introduce the parameter 
ot (see section 3.3) such that 


^'so. 


- ct*{> 


COS 


(3.4.6a) 


'H2SO4 


= - (1 -■ «Tt)>J 


COS 


(3.4.6b) 


These expressions for the flux of SO, and H2S0^ auto- 

maticallv satisfies the sulfur atom conservation relation. 

Eliminatina and 4 i„ from the oxygen atom 

CO 2 SO 2 "2“°4 

conservation relations, (3.4.6a) and (3.4.6b) we obtain 


= (5-2vx) - 2 


'*’cos 


TOS 


'^H20 (3.4.7) 


^COS 


For simplicity, consider that the hydrogen conservation 
relation is satisfied to first order by an upward flux of 
H 2 O and a downward flux of H 2 SO^. For this situation (3.4.7) 


becomes 


'*’co 

= (4 - a) + 2 

* 0 , I 

“ j 

‘^’cos 


'*’cos 


(3.4.8) 
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Consex'vation Relations 
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where and 


ative) and ^ 
case where X 


COS 


H2SO 


^CO/"COsi» 4, 


were assumed to be downward directed fluxes (nag- 

2 

an upward flux at the cloud tops. For the extreme 

» 1 (i.e. a * 0, » 0) at the cloud tops, 

2 ^2 

^C02 » 0, j'*'H2S0^/'^C0s[» 1 ,j’^H20/'^C0S )» 1 and 


j C02/^QQg| » 3 so that the over all chemical reaction below the 
cloud tops is characterized by 


4C0 + H2®°4 ^^°2 ^^2° (3.4.9) 


Th-^ chemical reaction was first suggested by Prinn (1973) for 
thermochemical formation of CO2, COS and H2O below the clo.ud layer. 


For the other extreme case >> 

“2 


Xso2 ^H2S0^ 


(i.e. I ^2/;^Qg 


"COS ■* "^S02 ' ^H2S04 ) at the cloud tops, j'^CO/^' 002 ) 1, |'02/?^q 

*■ so that the overall chemical reaction below the cloud layer 
is 


1 . 


CO + h02 -> CO2 (3.4.10) 

This is the effective reaction consistent with the modeling of 
Sze and McElroy (1975) . 


The real Venus photochemical model will produce results imme- 
diate to these two extreme cases. H2O molecules required in the 
formation of H2SO^ by reaction (3.4.3) may be either produced 
above the clouds or supplied by an upward flux from below the 
cloud. 
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CHAPTER IV 


RESULTS 


4 . 1 Pre-Sulfur Chemistry 

The atmospheric model used in this section is summarized 
in Table 1. The chemical model is presented in Table 3. This 
is similar to the model used by Sze and McElroy (1975) with the 
exception of several updated rate constant values. The eddy 

diffusion coefficient adopted here is derived from the equation 

* 5 2-1 

(3.1.2) with and K set equal to 2.5 x 10 cm sec and 

6 2 "X 

6 X 10® cm'^ sec” respectively. The boundary condition on the 0, 

. . -7 

mixing ratio at the cloud tops was set equal to 1 x 10 

The carbon and oxygen budget for the Venus atmosphere is 

1 1 

illustrated in Figure 3. Photolysis of CO 2 produces 1.5 x 10 

-2 -1 

molecules cm sec of carbon monoxide and an equal amount of 

atomic oxygen. Catalytic recombination above the cloud layer 

13 -2 -1 

returns 0.65 x 10 molecules cm sec of CO 2 . Carbon monoxide 

is catalytically converted into CO 2 through the two major paths 

simmarized by the net reactions (2.2.4) and (2.2.5). The excess 

13 -2 -1 

CO and 0 (0.85 X 10 cm sec ) will be transported below the 
cloud deck and recombined by unidentified chemistry, providing an 
upward flux of CO 2 to complete the carbon and oxygen cycles. 

The concentrations of key species calculated here differ 
from those of Sze and McElroy (1975) mainly because of revision 
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Corbon 










in rate constants. In regions near the cloud top, the OH con- 
centration, for instance, is about a factor cf 4 to 5 lower than 
that calculated previousli-. The HO 2 concentration is reduced by 
about a factor of 2, whereas the 0 and O 2 concentrations are 
enhanced by factors of about 3 and 2 respectively. The increase 
in 0 and O 2 abundances is caused by the reduction of odd hydrogen 
which results primarily from the larger rate constant for the 
reaction (2.2,9) as discussed in Section 2.2. The net efficiency 
for CO 2 catalytic recombination, however, is reduced by only a 
factor of 2 because the updated rate constant for the key reaction 

CO + on -> CO 2 + K (4.1.1) 

is about two times larger. 

4.2 Sulfur Chemistry 

The chemical model for the sulfur species (COS, S, SO, SO 2 , 

SO^, H2S0^) is given in Table 4. The photolysis of COS represents 

the sole source of sulfur considered in the present study. The 

source strength of sulfur is determined uniquely by the mixing 

ratio of COS at the cloud tops, for a given profile of the eddy 

diffusion coefficient. For the eddy diffusion coefficient 

specified in Section 4.1., photolysis of COS results in an in- 

12 -2 -1 

tegrated source for S of about 1 x 10 cm sec for a COS 

mixing ratio of lo”^ at the cloud top. Higher values for the 

eddy diffusion coefficient would result in a larger source of sulfur. 
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fhe sulfur atom released from COS photolysis will parti- 
cipate in a rapid reaction with O 2 

S + O 2 ^ SO + 0 , (4.2.1) 

forming SO and 0. This represents a significant source for 
atomic ©.xygen near the cloud top as illustrated earlier in 
Figure 2. 

Figure 4 shows the impact of sulfur chemistry on the key 
radial species 0, OH and HO^. The abundance of atomic oxygen 
below (38 km is significantly enhanced due to reaction (4.2.1). 

The increase in 0 conoantraticn results in a reduction of IIOj 
through the reaction 

0 + HOo •> OH + O 2 (4.2.2) 

It follows that formation H 2 O 2 through the reaction 

HO 2 + HO., .. H 2 O 2 + 0., (4.2.3) 

is suppressed significantly, resulting in a lesser production 
of OH through photolysis of H202« 

Figure 5 shows the distribution of the sulfur species. 
Carbonyl sulfide (COS) is the dominant sulfur species below 65 
km with SO 2 comparable to COS above this altitude. The abundances 
of the other sulfur species (S, SO, SO^) are considerably smaller. 

The formation of SO, is summarized in Figure 6. The dominant 
reaction is given by 
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Impact of Culfur Chemistry on Hio Nurnbor Pr'iisitv 
Profrlcs of O, Oil and IIO 2 . The mixing ratio„of 
O 2 at the cloud top was set equal to 4 x 10~®. 
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so + 0 + M -> SO 2 + M (4.2.4) 

which is an order of magnitude larger than reactions such as 

SO t O 3 > SO 2 + C>2 (4.2.5) 

SO + O 2 SO 2 + 0 (4.2.6) 

SO + OH SO 2 + H (4.2.7) 


Reaction (4.2.7) was thought by Prinn (1973, 1975) to be the 
dominant reaction in the formation of SO^. 

Figure 7 shows the production rate of SO^ , which is 
equivalent to the production rate of H,S 0 ^ molecules, as a 
function of altitude. The reaction 

SO 2 + 0 + M •> SO^ + M (4.2.3) 

provides the dominant source of SO^ and is about two orders of 
magnitude larger than the source derived from reaction 

SO 2 + HO 2 ->■ SO^ + OH (4.2.9) 

The reaction (4.2.9) was emphasized by Prinn (1973, 1975). 

It appears that the overall reaction scheme for K 2 S 0 ^ 
formation is represented by the reaction sequence (2.3.3) 
discussed in Section (2.3). This is essentially a sulfur and 
oxygen scheme in which odd hydrogen species are not directly 
involved. 
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One of the central issues in Venus aeronomy is understanding 
the fate of oxygen produced by CO2 photolysis. There are four 
paths for oxygen to follow. It may be catalytically recombined 
with CO to form CO2 as indicated by the net reaction schemes 
(2.2.4) and (2.2.5). Alternatively it may combine with sulfur 
forming SO2 and H2S0^ both of which would be transported down- 
ward to replenish the Venus cloud. The remaining possibility is 
the formation of O2 molecules which could also be tranported 
below the cloud deck. 


Figure 8 describes the apportionment of oxygen between these 

four paths as a function of COS mixing ratio at the cloud top. 

—8 

At low concentrations of COS ('v 10 (v/v) ) , catalytic recombination 

processes, ^^^t' ^ downward flux of O2/ , consume almost all 

the available oxygen provided by photolysis of CO2 • At larger 

concentrations of COS ("v 10~^ (v/v) ) , rapidly decreases with 

^2 

increasing COS. The fraction X^^^, however, decreases only 

slightly. The decrease in X_ is balanced by a corresponding 

increase in the downward flux of SO2 and ' <3escribed in 

Figure 8 by X__ and X„ -- respectively. 

" 2°°4 


For this particular calculation, the O2 flux is consumed 

_7 

above the cloud tops for a COS mxxing ratxo of about 4 x 10 
Ideally this optimizes the condition for formation of H2SO^. 
However, it may also lead to an SO2 abundance exceeding the 
observed upper limit. Computed values for X„,. , X^ , X_^ and 

Cat Urt 
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Impact of Sulfur on Oxygen Produced by CO 2 Photolysis. The dependence 
of the functions (see the text) and the colunm density of S02 on the 
abundance of COS is calculated for an O 2 mixing ratio of 4 x 10“^ at the 
cloud tops. 



HjSO^ as a function of COS mixing ratio are given in Table 8. 
Flux ratios of interest, evaluated at the cloud tops, are 
also presented as a function of the COS mixing ratio in 
Table 9. 


Given a downward flux of H2S0^ from the model calculation, 
the key question is how thick a cloud deck can be maintained 
under steady state conditions. Clearly the cloud abundance, N , 

w 

would be in part determined by the downward flux of H-5SO. , <{>« » 

and in part determined by the residence time, t, of the H2SO^ 
droplets in the cloud. The relationship between the three quantities 

Nc, gQ and t is given in Appendix A. The concentration of 

? ^ -5 

H2SQ^ in the cloud deck is assumed to have a mixing ratio of 10 

as suggested by Rossow and Sagan (1975) . We consider two diff- 
erent cloud thicknesses of 5 and 10 km and show in Figure 9 the 
relationship between __ and t . The residence time of the 
H2S0^ droplets is unknown, but is likely to lie betv/een 0.5 and 
3 years as indicated by the shaded area. 


For large x , about 3 years, a H2SO^ cloud as big as 10 km 

may be maintained by a COS mixing ratio as low as about 2 x lo”^. 

on the other hand, if the residence time is as small as 0.5 year 

or less, an exceeding large concentration of COS, approaching 
— 6 

10 , is required to maintain a 10 km thick cloud. However, there 

are additional constraints on the abundance of COS. The oxygen 

flux will be exhausted by sulfur species when the COS mixing ratio 

-7 

exceeds about 5 x 10 . In addition the calculated SO^ column 
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Fate of Oxyyen Atoms Produced by CO,^ Photolysi 


P 


<A 


T 


O i 

1 O 

fN 

C\ 

CN 

-5T 

CO 

[ >-1 

CO 

n 

C> 

O 


o 

o 

cN 

CN 

m 

i< 

• 

t 

• 

• 

• 


CI 

CO 

o 

rH 

o 

•T 

O 

o 

CO 


■T 

CN 

CO 

o 

o 

eg 

en 



• 

• 

• 

• 



CN 


CO 



CO 

o 


VO 

CO 

«T 

o 

X 

i in 


H 

O 

o 


f 

t 

• 

• 

• 

1 

XJ 

1 in 

o 

rH 

H 

o 


o 

r~ 

n 

IN 

r-i 

0 


01 

cn 

01 

01 


1 • 

• 

• 

• 

• 


o 

N 

w 

o 

u 

iw 


CO 

1 

i 

1 

r>- 


o 

1 

o 

1 

o 

1 

o 

1 

o 

rH 

H 

iH 

H 

rH 

t< 

« 


N 


H 

H 

01 

'3' 

in 


Ratio of Fluxes at the Cloud Toi>s 








w 

o 

u 


o 

•r 

m 

!N 

r> 

n 

in 


•T 

VO 

(n 


m 

in 

o 

• 

■ 

• 

• 

• 

o 

o 


r> 

fn 

Cl 

1 

r~i 

1 

1 

1 

1 


o 

•» 

m 

n! 

r» 

r> 

in 


■Sf 

VO 



n 

in 

o 

• 

• 

• 

• 

t 

68 

Ol 

n 

<N 

CN 


H 

CM 

CM 

Cl 


n 

n 

H 

o 

o 

(N 

o 

O 

o 

o 

• 

• 

• • 

t 

• 


w 

03 

00 

m 

■sr 

VO 

o 


VO 

00 

in 

in 

u 

iH 

n 

fO 

n 

n 

' 

« 

» 

« 

• 

« 


o 

cn 

CN 




in 


CO 

<N 

o 

in 

m 

m 

«5' 

o 

H 

o 

Cl 

VO 

VO 


u 



n 

m 

n 



• 

• 

• 

• 

t 



• • • • t 

I I I I I 



w 

o 

u 


nj 

O 

-> 


W 

o 

u 

■ 0 - 


00 

o 

in 


O' 

00 

o 

o 

01 

01 

in 

VO 

VO 

VO 

VO 

• 

1 

• 

1 

• 

1 

• 

1 

r 

«!r 

H 

CN 

tn 

n« 


iH 

r- 

00 

•4* 

01 

if) 

'T 

o 

H 

t 

• 

• 

• 

• 

01 

ro 

O 

o 

O 

"S' 

1 

1 

1 

) 

f 


N 


w 

o 

u 


4-1 


CO r* r« 

I I I I I 

o o o o o 

'■H H H H r-4 

X H X N X 

H H m 'T in 


Mixing Ratio of CO 



(j^) aujii ;uap|say ^OS^H 



abundance might exceed the observed upper limit for a COS mixing 

-7 

ratio larger than about 2 x 10 . The calculated SO 2 column 

abundance is however highly model dependent and critically 
relies upon the assumed eddy diffusion coefficient and the 
adopted chemistry scheme. Revision in rate constants or 
diffusion coefficients or omission of Key chemical reactions 
could modify the calculated distribution of sulfur species, 
particularly SO^ 

M 
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CHAPTER V 


CONCLUDING REMARKS 

We presented a model for Venus' atmosphere which com- 
plements and extends the earlier work by Prinn (1973/ 1975) 
and by Sze and McElroy (1975). We assumed that sulfur is 
released by photolysis of COS (Prinn, 1973) and that oxygen 
is supplied by photolysis of CO2. The supply rate for sulfur 
is determined by COS mixing ratio at cloud top and by dynami- 
cal processes. The oxygen supply rate, however, is influenced 
by a complex balance between odd hydrogen, sulfur and CO2 
photochemistry, in addition to transport processes. 


We confirm the importance of O2 as it controls the rate 
limiting step in the H2SO^ formation sequence. The presence 
of sulfur species may enhance the concentration of atomic 
oxygen above the cloud top due to the rapid reaction. 


S + O2 SO + 0 (5.1.1) 

with S released from photolysis of COS. Subsequent reactions 
with atomic oxygen, 

SO + 0 + M SO2 + M (5.1.2) 


and 


followed by 


SO2 + 0 + M ^ SO2 + M (5.1.3) 


SO3 + H2O H2S0^ (5.1.4) 

appear to be the dominant processes for H2S0^ production. 
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Reaction scheme (2.3.2) which involves odd hydrogen species 
plays a relatively minor role in the formation of H^SO^ 
molecules. The calculated SO2 abundance resulting from our 
present model appears to be quite high, but it could be 
lowered by using a faster rate constant for reaction (5.1.3) 
or by incorporating new reactions which transform SO2 to SO^ 
or H^SO^. Lower abundance of SO2 would imply a more efficient 
scheme for H2S0^ production. Future measurements on SO2 
abundance would place an extremely useful contraint on Venus' 
sulfur chemical model. 

For a given supply of 1 ^ 2^^ 4 n'oiscules, the abundance of 
sulfuric acid cloud should be controlled by complex dynamic 
processes. We oresented a steady state cloud model by intro- 
ducing the concepts of residence or characteristic time 

constants for Venus' cloud deck. Wa showed that a 10 V;m 

. -5 

cloud with a mean mixing ratio of H2SO^ 10 would require 
an extremely long residence time of about a few years. A 
short residence time ( < 1 year) would imply either a less 
abundant H2SO^ cloud deck or a calculated flux of H2SO^ that 
is too low. 

A significantly larger flux for H2SO^ than those cal- 
culated here would imply ?n alternate oxygen source for the 

oxidation of SO2 molecules. An example of such an oxidation 

1 3 

process is given by excited sulfur dioxide ( SO2, SO2 ) 
chemistry, resulting from photo-excitation processes. 
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+ hv ->• 


(5.1.5) 


SO,, 


'SO. 


and 


SO2 + hv -> '’SO2 


(5.1.6) 


and might extract oxygen directly from CO 2 viz 


CO 2 + SO 2 -i- CO + SO 3 (5.1.7) 


and 


CO 2 + SO 2 CO + SO 3 (5.1.8) 


Reaction scheme (5.1.5 - 5.1.8) would, in principle, 
bypass the limited supply of 02 * The potential role of 
excited sulfur chemistry in Venus’ cloud formation will be 
investigated in our future work. 
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APPENDIX A 


Relation Between COS Flux and Residence Time 
of H2S0^ Droplets for Different Cloud Thickness 


For a steady state cloud, the rate of supply of H2S0^ 
molecules to the cloud deck by downward flux, on ' must 
balance the rate of removal of H 2 SO^ cloud droplets \>rhich 
may be represented by 

N_ 


‘^H2S0^ 


(A.l) 


Where N^ is the column abundance of H2S0^ cloud and t is 
the residence time or the characteristic time constant for 
the cloud deck. 


The column density N^ depends upon the cloud thickness 
h, the average mixing ratio of H^SO. in the cloud, f„ „„ , 

** * rl fc) vj ^ 

and the mean background concentration of CO 2 , n^Q , in the 
cloud deck ; 


N, 


.-5 


Using the value 10 for f„ 

H2SO4 

19-3 £ ^ 

value of 1 X 10 cm for n^^^ , 

2 


^H2S04”c02^ 

, taken from Table 2, and a 


N, 


19 * 

1 X 10'"-^ h 


(A. 3) 
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*2 * 

where has units of cm and h is the cloud thickness 
expressed in km. 


The expression (A.l) may alternatively be written as a 
relation between t and the mixing ratio of COS at the cloud 
tops, flux of H 2 S 0 ^ at the cloud tops is 

given by 


'^ 250 ^ 


- (1-a) 


(A. 4) 


as discussed in Section 3.4. Furthermore the flux of COS 
at the cloud tops, determined uniquely by 

f^Os^^B^ for a given eddy diffusion profile. Using the eddy 

diffusion coefficient defined by (3.1.2) with and K set 

5 2 “I 6 2 —1 

equal to 2.5 x 10 cm sec and 6 x 10 cm sec , respect 

-ively, we have 


■^COS^^B^ 


1 Q 

1 X 10-^^ ^COS^^B^ 


where 

ing (A.l), 

is 


-2 -1 

(Zg) is expressed in units of cm sec . Combin- 
(A.3), (A. 4) and (A. 5), the desired relationship 

* 

h 




(1-a) 


(A. 6) 


where t is expressed in seconds. Assuming a representative 

* 

value for a of 0.5 and values of 5 km and 10 km for h , the 
results presented in Figure 9 were calculated from (A.l) , 
(A. 3) and (A. 6) . 


52 


References 


Anderson, R. C., Pipes, J. G., Broadfoot, A.L, and V^allace, L. , 

J. Atmos . Sci . , 874 (1969). 

Bader, L. W. , Ogryzlo, E. A., Nature , 201 , 491 (1964). 

Barker, E. S., Icarus, 25, 268 (1975). 

Baulch, D. C. , Drysdale, D. D. and Lloyd, A. C., High Temperature 
Reaction Data , 3 ^, 18 (1969) . 

Burrows, J. P., G. W. Harris and B. A. Thrush, Nature , 267 , 233 
(1977) . 

Chung, K. , Calvert, J. G. and Bottenheim, J. W. , Intern. J. Ghent . 
Kinet . , (1974) . 

dyne, M. A. A. and Watson, R. T., J. Chem . Soc . Fara . Trans . l., 
10, 2750 (1974). 

Clyne, M. A. A. and S. Down, J. Chem . Soc . Fara . Trans . II , 70 , 
253 (1974). 

Cruikshnznk, D. P., Commun . Lunar Planet Lab 6 , 199 (1967). 

Cruikshnzk, D. P. and Kuiper, G. P., Commun . Lunar Planet , Lab 6 , 
195 (1967). 


53 


t 


Davis, D. D., Klenun, R. B. and Pilling, M. , Intern . J. Chera . Kinetic s* , 
4, 367 U972). 

DeMore, W. , private communications (1976). 

Pair, R. W. and Thrush, B. A., Trans . Faraday Soc . , 65, 1557 
(1969) . 

Hack, W. W. Hoyerman and Wagner, H. G. , int . J. Chem . Kin . Symp . 

1, 329 (1975). 

Harker, A. B., Ph.D. Thesis, University of California, Berkeley 
(1972) . 

Hudson, R. D. , (ed. ) , "Chlorofluoromethanes and the Stratosphere." 

NASA Reference Publication 1010 (1977) . 

Johnston, H. S., Morris, E. D. , Jr. and Van den Bogaede, J. Chem . 

Soc . , 91, 7712 (1969) . 

Kaufman, F. , Elementary Gas Reactions . Annual Review of Physical 
Chemistry , 20 , 45 (1969). 

Kaufman, F., Ann . Geophys . , 20 , 106 (1964). 

Kumar, S. and Broadfoot, A. L., Geophys . Res . Lett . , 2, 357 (1975). 


54 


Marov, »M. Ya. , Avduevsky, V. S., Borodin, N. F., Ekonomov, A. P., 
Kerzhanovich, V. V., Lysov, V. P., Moshkin B. Ye., Rozhdestvensky, 
M. K. and Ryabov, 0. L., Icarus, 20, 407 (1973). 


KcCrimb, J. L. cuid F, Kaufman, J. Chem . Phy s . , 5" , 1-70 (1^72^ 

McElroy, M. B. and Donohue, T. M., Science , 177 , 986 (1972). 

McElroy, M. B., Sze, N. D. and Yung, Y. L. , J. Atmos . Sci . , 30, 
1437 (1973). 

Payne, W. A., Stief, L. J. and Davis, D. D., J. Amer . Chem . Soc 
95, 7614 (1973). 




Prinn, Ronald G., J. Atmos . Sci., 32, 1237 (1975). 

Prinn, Ronald G. , Science, 182 , 1132 (1973). 

Reeves, R. R. , Manella, G. and Hartsck, P., J. Chem . Phys . , 32 , 
632 (1960). 


Rossow, W.B. and Sagan, C. , J. Atmos. Sci., 1164 (1975). 


Schofield, K., J. Phys . Chem . , 2, 25 (1973). 

Sze, N. D. and McElroy, M. B., Planet Space Sci., 23, 763 (1975). 


Traub, W. A. and Carleton, N. P., J. Atmos . Sci., 1045 (1975) 


Vinogradov, A. P., Surkov, Vu. A., Andreichikov, B. M., Kalinkina, 
0. M. and Grechischeva, I. M. , Planetary Atmospheres , Reidel, 
Dordrecht, C. Sagan, T. C. Owen, H. J. Smith (eds.) p.3 (1971). 

Walker, R. , Ph.D. Thesis Unversity of London, Queen Mary Colleg 
(1972) . 

Young, L. D. G. , Icarus , 17 , 632 (1972). 

Zahniser, M. S. and Kaufman, F., j. Chem . Phys . , 66 , 3673 (1977). 


Zahniser, M. S., J. S. Chang and J. G. Anderson, Chem . Phys . Lett . 
27, 507 (1974). 


